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ABSTRACT 

The barred grand-design spiral M83 (NGC 5236) is one of the most studied galaxies given its 
proximity, orientation, and particular complexity. Nonetheless, many aspects of the central regions 
remain controversial conveying our limited understanding of the inner gas and stellar kinematics, and 
ultimately of the nucleus evolution. 

In this work, we present AO VLT-SINFONI data of its central ~ 235 x 140 pc with an unprece- 
dented spatial resolution of ~0.2arcsec, corresponding to ~4pc. We have focused our study on the 
distribution and kinematics of the stars and the ionised and molecular gas by studying in detail the 
Paa and Br7 emission, the H2 l-OS(l) line at 2.122 /im and the [Fell] line at 1.644 /im, together with 
the CO absorption bands at 2.293 //m and 2.323 /im. Our results reveal a complex situation where 
the gas and stellar kinematics are totally unrelated. Supernova explosions play an important role in 
shaping the gas kinematics, dominated by shocks and inflows at scales of tens of parsecs that make 
them unsuitable to derive general dynamical properties. 

We propose that the location of the nucleus of M83 is unlikely to be related to the off-centre 'optical 
nucleus'. The study of the stellar kinematics reveals that the optical nucleus is a gravitationally bound 
massive star cluster with Mdyn = (1-1 i 0.4) x 10^ M©, formed by a past starburst. The kinematic 
and photometric analysis of the cluster yield that the stellar content of the cluster is well described 
by an intermediate age population of logT(yr) = 8.0 ±0.4, with a mass of ^ (7.8 ±2.4) x 10^ M©. 

Subject headings: galaxies: individual (M83, NGC 5236) - galaxies:general - galaxies: starburst - galax- 
ies:nuclei - galaxies:kinematics and dynamics - galaxies: structure 



> 



O 
(N 



X 



1. INTRODUCTION 

M83 (NGC 5236) is a nearby {D = 4.6 Mpc, 
22pcarcsec-\ z = 0.0017 from the NASA/IPAC Ex- 
tragalactic Database, NED) barred grand-design spiral 
galaxy with a nuclear starburst. The galaxy has been 
object of intense study during the last decade, given 
the complexity of its central regions, its proximity and 
the fact that it is alm ost face-on, with an inclination of 
i = 24° ()Comtelll981f ). This makes it a good candidate on 
which to make use of high spatial resolution integral field 
spectroscopy (IFS) to study the controversial aspects of 
its innermost regions. 

The general morphology of the galaxy shows a pro- 
nounced bar and well-defined spiral arms where star for- 
mation is intense. On the other hand, the central regions 
of M83 in the infrared are rather complex. The J-K im- 
ages of the inner region show two non-concentric circum- 
nuclear dust rings, w hich are associate d with two inner 
Lindblad resonances (jElmegreen et al.1tl998,) . These two 
rings are connected by an inner bar, almost perpendic- 
ular to the main stellar bar. The general shape of the 
extended emission traces an arc between these two dust 
rings, where the star formation is concentrated. 

The location of the nucleus of M83 remains un- 
clear. iThatte et al.l |2000) first reported the existence 
of a 3.4arcsec (~ 75 pc) offset between the optical 
nucleus and the centre of symmetry of the bulge K- 
band isophotes. The centre of symmetry of these ex- 
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ternal isophotes is coi ncident wi th the dynamical cen- 
tre proposed by Sakam oto et al.l (j2004), based on two- 
dimensional CO spectrosco py, and confirm ed later by 
Rodri gues et al.l (|2Q09) andlKna^en et al.l (2010). Dif- 
ferent locations of hidden mass concentrations were pro- 
posed to host the supermassive black ho le of M83, mainly 
based on studies of t he gas kinematics (jMast et"ani2006L 
'Diaz e t ani2006L iRodrigues et aD[2009l and iKnapen et al.l 
2010). 

In this paper, we present new integral field VLT- 
SINFONI fei^enhauer et al.ll2003f ) spectroscopy in H±K 
bands, with an unprecedented spatial resolution, cov- 
ering the central ^ 235 x 140 pc of the galaxy. We 
study the stellar and gas kinematics of the inner parts 
and address some of the open questions regarding 
the hidden mass concentrations at of f -nuclear locations 
(|Thatte et al.l I2000L iMast et al.l I2006L iDiaz et al.l [2006f ) 
and the ori gin of the steep velocit y gradients in the gas 
kinematics (jRodrigues et al.l l2009f ) , revealing a complex 
scenario where supernovae play a key role in the kine- 
matics. 

2. OBSERVATIONS 

2.1. Observations and Data Reduction 

The M83 observations are divided into four different 
pointings, labeled as A, B, C and D in Fig. [TJ These 
pointings were chosen to cover the stellar nucleus of the 
galaxy, that is identified with the optical nucleus of M83, 
covered by pointing A; t he centre of symni etry of the 
bulge K-band isophotes ()Thatte et al.|[2QQQ[ ). that cor- 
responds to pointing B; the proposed location of a hid- 
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den mass by iMast et all (|2QQ6[ ) , by pointing C; and the 
putati ve massive black hole location given in iDiaz et aTl 
(|2QQ6l ) sampled by pointing D. The first three, A, B and 
C, were carried out between April and June 2009 in ser- 
vice mode, using the AO module fed by a LGS. The 
fourth was performed in July 2011, also in service mode. 
The data were taken in the H+K configuration, using a 
scale plate of 0.05 x 0.1 arcsecpixel"^ that yields a nomi- 
nal field of view (FoV) of -3.2" x 3.2" then enlarged by 
dithering. The wavelength range covered is from 1.45 /im 
to 2.46 /am with a spectral resolution of R'^1500. 

Although a total of four pointings were programmed, 
the main analysis in this work has been performed on 
three of them. A, B and C. The observations for point- 
ing D could not be completed and there was only 1 usable 
object frame. Taking into account the vast difference of 
quality in the data, we decided not to include the new 
data in the main analysis, but use it instead to support 
some of the results. The first three pointings cover an 
area of —8" x 6" around the nucleus of the galaxy, while 
the fourth pointing covers — 3" x 3". The footprints 
of the observed pointings are shown in Fig. [TJ The to- 
tal integration time was 3300 s for pointing A and 3600 s 
for each of pointings B and C, split into individual ex- 
posures of 300 s. In addition, four sky frames of 300 s 
were taken for every pointing every two on-source expo- 
sures to subtract the sky emission, following the pattern 
OOSOOSOO. In the same way, a total of five standard 
stars (Hip066957, Hip069230, Hip070506, Hip071136 and 
Hip098641) were observed in order to perform telluric 
and fiux calibration. The fourth pointing is a single on- 
source exposure of 300 s with a matching sky frame. The 
standard star used for the calibration was Hip001115. 

The reduction of the data was performed using the 
standard ESO pipeline. An individual cube was built 
from each frame, from which the background sky emis- 
sion w as subtracted using the method outlined in Davies 
()2007[ ). We performed the telluric and fiux calibration 
on each cube individually to improve the results. Taking 
into account the relative shifts of the jitter pattern, we 
then combined the data to create a single cube for each 
pointing. Finally, these were combined to build a single 
mosaic. 

The telluric and fiux calibration were performed in 
two steps. First, each individual star was normalised 
to the continuum level, using a blackbody profile at 
the Teff listed in t he Tycho-2 Spectral Type Catalog 
(j Wright et allHool ). To remove the absorption features 
in the spectra of the stars, we used a solar template, con- 
volved and binned to match the resolution of the SIN- 
FONI data. The result is a "sensitivity function" that 
takes into account the atmospheric transmission. Sec- 
ondly, we used the H and K magnit udes of the stars 
from the 2MASS catalog (.Skrutskie et al.ll2006i) to con- 
vert our spectra from counts to physical units. We made 
us e of the response cu rves of 2MASS filters, as defined 
in iCohen et al.l ()2003l ). to obtain the values in counts 
of our spectra at their effective wavelengths. Using the 
above mentioned magnitudes, we obtained two conver- 
sion factors for each star (one for H-band and one for 
K-band). These factors were almost identical for both 
bands in every star, which verifies the calibration. We 
adopted the mean value to scale our curves. The fiux- 
calibrated cubes were obtained by dividing the individual 
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Figure 1. Central region of M83: HST NICMOS F222M image 
with the fields covered by the four SINFONI pointings superim- 
posed. The location of the optical nucleus and the arc of star 
formation mentioned in the text are also shown. Although the in- 
ner arc is not visible in the continuum map, it is clearly visible in 
the western part of the FoV in the flux panels in Fig. 3] The total 
coverage of our SINFONI observations is ~ 235 x 140 pc. 

cubes by the "sensitivity function", to correct from the 
atmospheric transmission, and by multiplying them by 
the conversion factor. 

2.2. Gas and Stellar Kinematics 

The gas kinematics were extracted by fitting a gaus- 
sian profile to the most releva nt emission l i nes, usin^ 
the code LINEFIT described inlDavies et al.l (|2Qllf ) (see 
also iForster-Schreiber et al.l l2QQ9f ). During the extrac- 
tion, the OH sky line at 2.18 /im is used to remove the 
instrumental broadening, measured to be 13 A FWHM. 

To extract the stellar kinematics, we focused on 
the two most prominent CO bands, CO (2 - 0) at 
2.293 /im and CO (3 - 1) at 2.323 /im, and used the 
Penalized Pixel-Fitting (pPX F) software developed by 
iCappellari fc EmsellemI (|2QQ4l ) to fit a library of stellar 
templates to our data (see Fig. [2j). We made us e of the 
near-I R library of spectral templates from Winge et al.l 
(120091) . which covers the wavelength range of 2.15 /im - 
2.43/im with a spectral resolution of R ^ 5600 and sam- 
pled at lApixel"^. The library contains a total of 23 
late-type stars, from F7III to M3III, and was previously 
convolved to our SINFONI resolution. 

The uncertainty in the gas and stellar kinematics is 
highly dependant on the S/N and how well resolved the 
line is. For the gas kinematics, the uncertainties could 
range from ^ 1% of the resolution element in those re- 
gions with high S/N up to more than ~10% in the re- 
gions with poorer S/N. The kinematic precision achiev- 
able with the stellar absorption features also depends on 
the S/N, although it is typically lower than the precision 
in the gas kinematics mainly due to uncertainties in the 
template matching uncertainties. As shown in Fig.[2l the 
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Figure 2. Integrated spectra of apertures A, C, E and H, rep- 
resented in Fig. and the result from fitting template spectra to 
the stellar absorption features. The normalised spectra are shown 
in black, and the results from the fitting are plotted as a thick blue 
line. The red dotted lines show the residuals from subtracting the 
fit. 

quality of the template fitting indicates that the preci- 
sion we achieved in the stellar kinematics is high, less 
than '^40kms~^ in those regions with high S/N, and 
that offsets of >30kms~^ would be clearly visible in the 
residuals of the fitting. 

In order to compare the velocity fields of the different 
phases of the gas and the stellar component, we have es- 
tablished a reference value of cz = 589.6 km s~^ for the 
velocity that has been used as a zero-point for the veloc- 
ity maps. This reference value for the velocity has been 
chosen as the mean value of the stellar velocity in a small 
aperture of 5 spaxel radius ce ntred in the dynamical c en- 
tre of the galaxy proposed by iSakamoto et al.l ()2QQ4l ) . 

2.3. Voronoi Binning 

Before the extraction of the kinematics, the data were 
binned to achieve a minimum value of S/N on the whole 
fie ld. We used the Voronoi bin ning method implemented 
bvlCappellari fc CopinI (|2QQ3f ) to maintain the maximum 
spatial resolution of our maps while constraining the min- 
imum S/N ratio. 

The Voronoi tessellation employs bins of approximate 
circular shape to divide the space, described in terms of 
a set of points called generators. Every bin encloses all 
the points that are closer to its generator than any other 
generator. The algorithm finds an initial set of genera- 
tors by selecting the spaxel with highest S/N ratio and 
accreting spaxels to that bin, until the required thresh- 
old is reached. Then, it moves downwards to lower S/N 
spaxels until all the points are assigned to a proper bin. 
This set of generators is then refined to satisfy both topo- 
logical and morphological criteria and to ensure that the 



Table 1 

S/N thresholds used for 
the Voronoi binning 



Feature S/N 



Paa 50 

Br7 45 

H2 l-OS(l) 15 

[Fell] 25 

Stellar 50 



scatter of the S/N among all the bins is reduced to a min- 
imum (see Cappellari & Copin 2003 for further details). 

This binning of the data does not affect those spax- 
els with high S/N ratio, preserving the original spatial 
resolution of these regions. The maps obtained for each 
line and the continuum are then binned independently, 
since their flux distributions are totally different. We 
have therefore defined different S/N thresholds for each 
Hue, in order to obtain appropriately sampled maps. The 
S/N cutoff used for each Hue and for the stellar contin- 
uum are shown in Table [TJ We selected these values to 
achieve roughly the same number of bins in each map. 

2.4. PSF Determination 

The measurement of the PSF has been done by com- 
paring the SINFONI data to higher resolution HST NIC- 
MOS data. As discussed in Davies (2008), we can use a 
higher resolution image with a well known PSF to esti- 
mate the PSF of a lower resolution image. After resam- 
pling the NICMOS image to our SINFONI pixel scale, 
the aim is to find a broadening function, 5, that sat- 
isfies IsiNFONi = Inicmos ^ B. Since the NICMOS 
PSF is weh known, we can estimate the SINFONI PSF 
as PSFsiNFONi = PSFnicmos ^ B. The shape of the 
resulting PSF is dominated by the broadening function. 

We performed independent fittings of the broadening 
function for each pointing and each band. The NIC- 
MOS images were obtained with the NIC2 camera, using 
the F160W and F222M filters for H and K band respec- 
tively, with a pixel scale of 0.075 arcsec pixel" ^. After 
trying different models for the broadening function 
we find that it is better described as a symmetric double 
Gaussian with a narrow component of ~1 pixel FWHM 
and another wide but fainter component of ^4-6 pixel 
FWHM that takes into account the seeing-limited halo. 
The different values of the FWHM of the resulting PSF 
are shown in Table O 

We note that the resolution of our data is limited by 
the pixel scale chosen for the observations rather than 
the LGS-AO performance since, in order to cover a wider 
FoV (~3arcsec), we chose the OVOSxOVl pixel scale from 
the three configurations available for SINFONI. 

3. OVERVIEW OF DATA 

The inner ~190 x 130 pc of M83 are covered by point- 
ing A centered on the optical nucleus of the galaxy, point- 
ing B on the photometric centre, and pointi ng; C on the 
off-nu clear black hole location proposed by iMast etld] 
(|2006h . 

The wide spectral coverage of the H+K band config- 
uration of SINFONI allows us to study in detail a large 
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Table 2 

FWHM of the PSF for each pointing. The values are measured after convolving 
the broadening function, B, obtained from the fitting with the PSF of NICMOS 

images. 



Pointing 


H-Band FWHM 




K-Band FWHM 




(pixel) 


n 


(pc) 


(pixel) 


n 


(pc) 


A 


3.66 


0.18 


4.03 


4.05 


0.20 


4.45 


B 


3.76 


0.19 


4.14 


4.18 


0.21 


4.60 


C 


2.95 


0.15 


3.24 


4.09 


0.20 


4.50 



number of emission lines and stellar absorption features 
(see Fig.[3|). In order to achieve a good level of S/N in the 
whole FoV, we focussed our study of the gas kinematics 
on the brightest emission lines, i.e. Br7 2.166 /im for the 
ionized gas, the ro- vibrational transition H2 l-OS(l) at 
2.122 /im for the warm molecular gas and the [Fell] line 
at 1.644 /im. We have extracted surface brightness, ve- 
locity dispersion and velocity maps of these three lines, 
represented in Fig. IH that allow us to study different 
phases of the interstellar medium. 

Although the Paa line is the brightest in the wave- 
length range covering the H and K bands, it lies at a 
wavelength where the atmospheric transmission is very 
low. Here, the strong variability of the sky absorption 
makes it very difficult to perform a good correction of 
the transmission, which translates into an increase of 
the noise compromising the results of the kinematics ex- 
traction. However, given the brightness of the Paa line, 
we were able to obtain a map of the emission that can 
be compared with HST/NICMOS archive images of the 
same region. 

As noted above, we make use of the H2 l-OS(l) line to 
trace the warm molecular gas in the whole FoV. How- 
ever, the detection of additional H2 transitions allows us 
to study in more detail, in Sec. HI the excitation mech- 
anisms of the H2 in the inner regions of M83, and dis- 
tinguish collisional excitation in shocks from radiative 
fluorescence. To improve the S/N ratio of the weaker 
transitions, we have integrated the signal of all the spax- 
els from the inner star-forming arc (those above a certain 
flux level) and those from the optical nucleus (see Fig.[T]). 
This allows us to measure the fluxes of six different tran- 
sitions with good level of confidence. The different lines 
measured are listed in Table [3l 

In addition to the numerous emission lines available, 
we detect a variety of stellar features in our spectra, 
most notably the Nal, Cal and CO absorption bands. 
As mentioned previously, we have focussed on the first 
two CO bands, CO (2 - 0) at 2.293 /im and CO (3 - 1) 
at 2.323 /im to study the stellar kinematics (Fig. [Sj). 

3.1. Morphology and Kinematics of the Gas 

As shown in Fig. 21 the morphology and, to some ex- 
tent, the kinematics of the distinct phases of the gas are 
rather different. The Br7 emission is mainly associated 
with the inner star-forming arc towards the west side 
of the FoV. The global velocity gradient of ~60kms~^ 
from northwest to southeast is consistent with an in- 
flow of gas alon^ the spi ral ar ms and through th e in- 
ner bar ([E lmegreen et_aLl [19981 Erosthwaite etlD[2QQ2|. 
iFathi et al..,2QQ8.) to the photometric centre of the galaxy. 



Superimposed on this, there is a ring that is on and 
around aperture A in Fig. HJ One of the most remark- 
able properties of this ring-like feature is that it shows no 
velocity gradient and has a low velocity dispersion that, 
together, argue against it being a dynamical structure. 
We discuss this feature in detail in Sec. [5l 

The H2 emission is mainly associated with the inner 
arc along the western part of the FoV and with the op- 
tical nucleus. The kinematics are very similar to those 
exhibited by the ionized gas. Around region A in Fig. 21 
the velocity is the same as that of the Br7 emission, 
where no dominant velocity gradient is observed. More- 
over, the emission of the molecular gas resembles the 
observed ring-like feature of the ionized gas. However, 
the strong emission of the optical nucleus allowed us to 
trace the velocity field across this region in more detail. 
The velocity gradient measured from northeast to south- 
west is ~200kms~^ in ~45pc, significantly steeper than 
the gradient measured from the stellar kinematics (see 
FigEJ. The velocity dispersion shows a similar picture 
as the Br7 emission, where the low values measured along 
the inner arc suggest that the gas is confined to a thin 
plane, presumably a disk supported by rotation. 

The [Fell] emission is highly extended along the inner 
arc, showing various knots of strong emission. Two of the 
brightest spots are located in the outskirts of the optical 
nucleus, labeled as E and F in Fig. IH at radial distances 
of ~ 30 pc and ~ 32 pc respectively. These two knots also 
exhibit a high velocity dispersion (~ 150kms~^ in region 
E and ~ 115kms~^ in region F, both taking into account 
the spatial resolution of our data) and are probably trac- 
ing individual super novae. We return to this issue in 
Sec. [5l Although the velocity field of the [Fell] is very 
similar to those traced by the Br7 and H2 emission, the 
velocity dispersion is systematically higher along the in- 
ner star- forming arc. As we discuss in Sec. [51 this higher 
velocity dispersion may be a sign of recent supernovae 
explosions, and would set a constraint on the age of the 
stellar populations along the arc. 

3.2. Stellar Component 

The stellar continuum derived from the first two CO 
absorption bandheads is mainly concentrated in the op- 
tical nucleus and, unsurprisingly, shows a similar mor- 
phology to that of the K-band image. On the other hand, 
the stellar kinematics show a completely different picture 
to that traced by the gas. Whereas the gas kinematics 
appear to be completely dominated by shocks and out- 
flows at small scales, the stars show a smooth velocity 
gradient from northeast to southwest, typical of a rotat- 
ing disk. Superimposed on this, the continuum emission 
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Figure 3. Normalised H+K spectra of the apertures A, C, E and H (see Fig. 3]for reference). Aperture A is located at the maximum of 
the Br7 emission, aperture C covers the centre of the optical nucleus, aperture E corresponds to one of the bright spots of [Fell] emission 
next to the optical nucleus and aperture H is located at the position of one of the SNR listed in Dopita et al.. (,2010 V The wavelengths of a 
number of lines and features are identified, and those studied in this paper are identified in blue. These spectra clearly illustrate the wide 
variety of excitation conditions that are occurring in the inner regions of M83. 



from the optical nucleus is dominated by a coherent in- 
ternal rotation, as highlighted in Fig. [6l The amplitude 
of the projected velocity field, measured peak to peak, is 
~ 32kms~^ within ~ 24 pc. The de-projected rotation 
velocities for an adopted inclination ofi = 24° would be a 
factor 2.5 higher. Although the uncertainties in the stel- 
lar kinematics are comparable to this value (in an spaxel 
basis). Fig. [6] shows a very clear jump in the projected 
velocity (magenta line) across the optical nucleus. 

As shown in Fig. [5l we detect two bright spots of con- 
tinuum emission (apertures H and I) in the outskirts of 
the optical nucleus, with a substantially high equivalent 
width. Given the spatial resolution of our data, these 
bright sources in the stellar continuum could be identi- 
fied as individual stars - late type giants or supergiants 
- in a post main sequence phase. This is consistent with 
the main scenario proposed for the optical nucleus, dis- 
cussed in Sec. [6l where the UV photons of a population 
of non- ionising stars would excite the molecular gas, ex- 
plaining the overpopulation of the Ji = 3 levels shown in 
Fig. [H Such a population of stars would be consistent 
with the age of the cluster derived from the CO absorp- 
tion bands, as discussed in SecEl 

3.3. Pointing D 

Because of the different integration time for the SIN- 
FONI data from pointing D, we have not included them 



in the general analysis of the emission and kinematics 
of the nuclear regions of M83. However, we can draw 
some con clusions about the hidden mass location pro- 
posed by ' Diaz et aTl ()2QQ6f ). As shown in Fig. [71 we 
have extracted a K-band continuum image, H-K colour 
map and Br7 surface brightness and velocity maps of 
the ~ 65 X 65 pc region. The H-K values obtained lie 
within the range ^ 0.2 — 0.7, i n good agreement with 
those derived by I Wiklind et all (j2004f ) using NICMOS 
F160W and F222W images. The dust lane that crosses 
from north to south in the central region of M83 is just 
at the east of our field of view, although the extinction 
gradient is clearly visible towards the western part of the 
map. The Br7 maps show the north-west end of the ring 
feature seen in Fig. IH and the smooth velocity field that 
continues beyond pointing C. The velocity gradient does 
not show any evidence of a hidden mass in this position, 
but simply refiects the difference in velocity between two 
regions, the ring-like structure observed in pointing C 
and the bright lane of Br7 emission at the west of point- 
ing D (see Fig. [7]). 

4. WARM MOLECULAR GAS: H2 TRANSITIONS 

There are two very distinct areas of strong H2 emission, 
one associated with the optical nucleus and the other 
with emission from the inner arc at the western part of 
the FoV (see Fig.|4]). We have extracted integrated spec- 
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Figure 4. Flux and kinematics maps of the main emission lines. From top to bottom, the Brj , H2 l-OS(l) and [Fell] maps, and 
from left to right, flux, velocity dispersion and velocity fleld. The boxes indicate the apertures used to extract spectra, some of which 
are shown in Fig. O These also act as reference points with respect to the discussion in the text. Flux maps are scaled with a factor 
3xl0~^^ ergs~^ cm"^. 



7 



Continuum 




^ 01 



- 137 



95 



RA (orcsec) 




RA (orcsec) 



14.25 



11.50 



8.75 



6.00 
70 



15 



-39 



-94 



-150 



Figure 5. Maps of the stellar continuum flux (top left), equivalent width (top right), velocity dispersion (bottom left) and velocity fleld 
(bottom right) for the nuclear region of M83. These were obtained by fltting the flrst two CO absorption bands, CO (2 - 0) at 2.293 /xm 
and CO (3 - 1) at 2.323 /im. The boxes indicate the same apertures as in Fig. |4] Flux map is scaled with a factor 3xl0~^^ ergs"-*^ cm~^. 
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Figure 6. Flux (white), velocity dispersion (green) and velocity 
(magenta) profiles along a pseudo-slit traced to include the optical 
nucleus along the direction of maximum variation of the velocity 
field. The location of the photometric centre of Thatte et al. (2000) 
is marked as a red dot. The pseudo-slit is plotted in black over the 
stellar continuum map for reference. 



tra from both the optical nucleus and the arc, including 
all the spaxels from each region above a given flux thresh- 
old. This threshold has been chosen to be 15% of the 
brightness of the H2 peak, and allows us to reject most 
of the weakest spaxels that contribute mostly to increase 
the noise. We have measured the fluxes of the differ- 
ent transitions by fitting a Gaussian profile, and derived 
the uncertainties using a Monte Carlo technique. The 
method consists of measuring the noise in the spectra as 
the rms of the residuals after subtraction of the Gaussian 
profile. Taking into account this estimation of the noise, 
we construct a total of TV = 1000 simulations of our spec- 
tra where the lines are again fitted. The uncertainty of 
our measurements is defined as the standard deviation 
of the fluxes for each line. The values obtained for the 
fluxes of the different transitions and their uncertainties 
are Hsted in Table [H 

Using these fluxes, we can calculate the various popu- 
lations of the upper H2 levels associated with each tran- 
sition. As shown in the population diagrams in Fig. [H 
we found that the emission on both regions has an im- 
portant contribution from non-thermal processes. We 
have compared the population of the different levels of 
both regions with three of t he photon- dominate d region 
(PDR) models discussed in iDavies et al.l (|2003 ^. where 
the excitation is dominated by far-UV photons. In Fig. [3 
model 1 and 2 consist of moderate and high density 
gas, nH = lO^cm"^ and nn = 10^ cm~^ respectively, 
where the main excitation process is UV fluorescence, 
whereas model 3 describes a fully thermalised region with 
nn 10"^ cm-3 and T = 2 x 10^ K. 

The emission from the inner arc shows that the low- 
est (z^ = 1) transitions are thermalized while the higher 
(z/ = 2) transitions are slightly overpopulated. This is a 
clear sign of fluorescent excitation mechanisms, that tend 
to excite the highest levels. The values are consistent 
with model 2 of a dense PDR in which the u = 1 levels 
are thermalised and the u = 2 levels are overpopulated 
by fluorescence. Intriguingly, the optical nucleus shows 
a stronger contribution of radiative processes, most no- 
tably in terms of an overpopulation of the = 5 levels 
(equivalently stronger S(3) lines). 

Although these results do not allow us to quantify the 
contribution from the different mechanisms, it is clear 
that the dominant processes are rather different in both 



Table 3 

Integrated fiuxes of the H2 lines for the optical 
nucleus and the inner arc. 



Line 


A 


Flux 






(/im) 


(10-15 ergs- 


icm-2) 






Optical Nucleus 


Inner Arc 


1-0S(3) 


1.958 


4.24 ±0.73 


1.19 ±0.20 


1-0S(2) 


2.034 


0.50 ±0.27 


0.43 ±0.06 


2-lS(3) 


2.073 


1.30 ±0.19 


0.23 ±0.05 


l-OS(l) 


2.122 


1.85 ±0.26 


1.24 ±0.15 


l-OS(O) 


2.223 


0.73 ±0.14 


0.44 ±0.05 


2-lS(l) 


2.248 


0.55 ±0.19 


0.34 ±0.07 



regions: while the inner arc seems to be dominated by 
thermal processes compatible with episodes of recent su- 
pernova activity, the contribution from radiative pro- 
cesses in the optical nucleus associated to fluorescent ex- 
citation mechanisms is highly significant. 

5. EMISSION AND KINEMATICS OF THE GAS: EVIDENCE 
FOR SUPERNOVAE 

As shown in Fig. HI the kinematics of the gas is totally 
unrelated to the stellar kinematics (Fig. [5]). There is 
no clear evidence of a single uniform velocity gradient, 
as seen for the stars, in any of the emission line maps. 
Together with the dispersion, this suggest that the gas 
kinematics are at small scales dominated by shocks and 
flows. In this section we argue that these characteristics 
are related to the presence of super novae. 

The nature of the Br7 ring-like feature noted in Sec. 13. II 
could be explained in terms of a light echo from a recent 
type II supernova explosion. This would be consistent 
with the low velocity dispersion and the presence of the 
complete Brackett series in the spectrum (see the spec- 
trum of aperture A in Fig. [3]). The projected radius of 
the ring is ~ 23 pc, which means that, if its origin was 
a supernova event, the explosion would have occurred 
~ 75 yr ago. If it is indeed a light echo, one might expect 
to see changes on a timescale of 10 years, for example 
between the NICMOS data from 1998 and our data from 
2009. We have compared the Paa luminosity profiles of 
the ring measured by NICMOS with our data to confirm 
a possible evolution with time. We first matched the 
PSFs, and then extracted horizontal and vertical profiles 
centred in the centre of symmetry of the ring feature in 
both images. We applied a single scaling and normalisa- 
tion to the profiles, derived to match the background of 
the emission far from the Paa ring. The two resulting 
profiles of the Paa emission are shown in Fig. [9l Consid- 
ering the expansion rate of a light echo, the difference in 
size of the ring between the two datasets is expected to be 
of ~3.4pc, which is less than a resolution element in our 
coarse sampling (see Table [2j). However, there is a signif- 
icant relative decrease of ~6% in brightness between the 
two epochs that cannot be explain as a PSF or AO effect. 
This decrease in the emission supports the hypothesis of 
a transient event like a supernova explosion. 

In Sec. 13.11 we also pointed out two bright [Fell] spots 
located at the outskirts of the optical nucleus, labeled as 
E and F in FigHl which are associated with high veloc- 
ity dispersion and are probably tracing individual super- 
novae. The spectrum of one of these sources (E) is shown 
in Fig. [3l The region labeled as D in Fig. |4] corresponds 
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Figure 7. Maps of the K-band continuum emission (left), Br7 surface brightness (centre) and velocity (right) from SINFONI data at 
pointing D. The white contours in the left panel correspond to the H-K colour obtained from the SINFONI data; the contours in the centre 
and right panels show the K-band continuum emission. The values of the Br7 velocity are in the same colour scale as those shown in Fig.^] 
for the other three pointings. 
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Figure 8. H2 excitation diagrams relative to l-OS(l) for the optical nucleus (left panel) and the inner star-forming arc (righ t panel). The 
lines from which the population levels are derived, are indicated. Overplotted are three of the five FDR models discussed in Davies et al.l 
{2003). Models 1 and 2 consist of moderate and high density gas respectively where the main excitation process is UV fiuorescence, while 
model 3 describes a fully thermalised region. As mentioned in the text, the error bars are obtained by a bootstrap method of = 1000 
simulations of the spectra. The diagram of the optical nucleus shows a strong overpopulation of the Jn = 5 transitions that points to 
non-thermalized excitation mechanisms like radiative fiuorescence, whereas the populations of the inner arc are consistent with a dense 
FDR where the u = 1 levels are thermalized by collisions while the u = 2 overpopulation is characteristic of fiuorescent excitation. 



to the source M83-SNR-N-01 identified in 'Do pita et al l 
(|2Q1Q| ) in their Table 3, and shows some characteristics 
expected for a recent supernova, i.e. strong [Fell] and 
H2 emission. 

There are also a few other supernova remnant (SNR) 
candidates that lie within our FoV, but none of them are 
obviously detected. The reason why there is no sign of 
most of these sources is because the selection criteria of 
the SNR candidates in lDopita et al.l ()2Q1Q| ) is in terms of 
the [Oil] emission at 3727 A and 3729 A. To achieve high 
[Oil] /H a ratios a radiative shock of ^300 -500 kms~ ^ is 
needed (|Dopita Sutherland! 119951 . Il996h . so it is the 



ionized pre-shock that actually emits in [Oil]. Therefore, 
the temperature and the shock speed ar e too high to ex- 
pect [Fell] or H2 emission fsee lBurton et all 119901 ). We 
have also as sessed the list of individual X-ray sources of 
iSoria fc Wul ([20b2) but only three of them lie within our 
FoV (sources 37, 40 and 43). Source 43 corresponds to 
the optical nucleus, source 37 is shifted by ^0.7^' to the 
west and by ~0.5'' to the south of our aperture A and 
source 40 is shifted ^0.5^' to the west and ^0.6" to the 
south of aperture B. Taking into account the uncertain- 
ties in the position of the X-ray sources, source 37 could 
be tentatively associated with the ring of Br7 emission. 
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As proposed by iRavmondl (|2QQl l), the width of the 
[Fell] and H2 Hnes is a good estimation for the velocity of 
the shock speed in SNRs. The [Fell] line is clearly broad- 
ened with an intrinsic width of ~100-180kms~^. In con- 
trast, the H2 line width is no more than ~40kms~^, in- 
dicating it is barely resolved. Thes e values are in good 
agreement with the values derived in lBurton et al.l ([l990) 
for fast J shocks (i.e. "jump" shocks, exhibiting discon- 
tinuous transition; in contrast to "continuous" shocks 
showing continuous transition in velocity, density and 
temperature) where the [Fell] emission is expected to be 
stronger, and slow J or fast C shocks for the molecular hy- 
drogen emission. Given the size of these regions of ~4 pc 
or less and the velocity dispersion of ~100-180kms~^, 
and assuming a constant expansion rate, we obtain an 
upper limit to their age of ~ 4 x 10^ yr which supports 
the argument that they are recent events. 

The global velocity gradient in the [Fell] emission of 
^60kms~^ along the arc is very similar to that observed 
for the Br7 line. However, the extended [Fell] emission 
throughout the inner arc and the high velocity dispersion 
are rather different compared to the Br7 and H2. Both 
are consistent with a scenario where the most massive 
stars have already exploded as supernovae. These explo- 
sions would blow up the gas of a thin quiescent disk in 
the perpendicular direction. 

This would provide a natural explanation for the sys- 
tematically higher velocity dispersion observed for the 
[Fell] emission along the inner arc: the shocks of the 
supernovae would be fast enough to dissociate the H2 
molecules and enhance the gas-phase Fe abundance as 
well as to generate singly ionised Fe, but not to fully 
ionise the H. 

A scenario in which many stars have recently exploded 
as supernovae, or may soon do so, is supported by th e 
strong radio continuum found by iSaikia et al.l () 19941) 
and the diff use X-ray emission detected in the arc (see 
iSoria fc Wul [2002). These authors also found a high 
abundance of Ne, Mg, Si and S with respect to Fe. This 
suggest that the interstellar medium would have been 
enriched due to type-II supernova explosions in the re- 
cent past. This scenario would also provide an origin for 
the radiative fluorescence excitation of the H2 molecules. 
It can be explained in terms of a population of stars of 
~5-8 M0 that have not exploded yet as supernovae (and 
may not do so), but emit e nough UV radiation to excite 
the H2 (^Puxlev et al.lll990f ). 

This scenario would also be in good agreement with 
the ages of optica lly-selected star clusters inferred in 
iHarris et al.l (|200Th . Based on HST WFPC2 observa- 
tions, these authors studied the star formation history of 
the southern star-forming arc and found that more than 
'^75 % of the more massive clusters (M > 2x 10^ Mq) have 
ages less than 10 Myr. They found a sharp cutoff in the 
age distribution of the star clusters, and proposed a re- 
cent burst of star formation activity that began < 10 Myr 
ago, suggesting an outward propagation across the arc. 

We can make a rough estimate of the age of the 
youngest stellar population along the inner star-forming 
arc under the assumption that the [Fell] emission origi- 
nates in supernovae. T he STARBURS T99 stellar popu- 
lation synthesis models (jLeitherer et al . 1999, SB99 here- 
after) predict that this emission would reach its maxi- 
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Figure 9. Flux profiles in arbitrary units of the Pao; ring-like 
emission from NICMOS and SINFONI. Distances are measured 
from the centre of the ring. Both images were previously PSF- 
matched and normalised to the background of the emission before 
comparing the data. There is a significative relative decrease of 
~6% in fiux between the two epochs. 

mum at ~10Myr (see Fig. [TQ|) . On the other hand, the 
Br7 emission also appears to be dominated by a super- 
nova event. It shows little ionized gas emission from the 
star themselves and it is probably part of the tail end of 
the supernovae for the same burst of star formation that 
will be discussed in Sec. [6l This implies that there are 
few OB stars left and that the star formation episode was 
consequently a short burst. Compared to the arc, the 
optical nucleus shows more H2 emission and less [Fell] 
while still having very little Br7 emission. This is also 
suggestive of a short burst of star formation. Both re- 
gions therefore appear to have experienced an episode of 
star formation around >10Myr ago. More detailed as- 
sessments of the age of the star formation in the optical 
nucleus, based on dynamics as well as spectral features, 
are presented in Sec. [6l 

6. OPTICAL NUCLEUS: AN EVOLVED MASSIVE 
OFF-NUCLEAR STAR CLUSTER 

In order to shed some light on what one should con- 
sider to be the nucleus of M83, we have studied in more 
detail the kinematic and photometric properties of the 
optical nucleus. We made use of NIC2 F222M K-band 
continuum image to fit the brightness profile of the clus- 
ter. We adopted a model of a symmetric Sersic profile 
convolved with the PSF to fit the core of the nucleus, 
combined with an asymmetric Gaussian to take care of 
the extended emission. The fits reveal a strongly peaked 
nucleus of i^eff = 2.97 ± 0.15 pc (with a Sersic index of 
n = 2.7), clearly unresolved in the SINFONI data, which 
sits on a more extended emission of size ~50pc. Given 
that the core of the cluster appears point-like and is not 
resolved in the AO data, we use this value as an upper 
limit for the effective radius of the cluster. 
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Figure 10. Evolution of the Btj equivalent width (left) and [Fell] flux (centre) from SB99 models, and C02.22-2.37/^m index (right) from 
M05 models. The dotted l ine in t h e left panel shows the upper limit of the Br7 equivalent width at the optical nucleus. We have assumed 
the theoretical ratio from [Colinal ( |1993| ) to convert the [Fell] 1.26 /xm flux predictions to [Fell] 1.64 /xm fluxes. The upper limit for the 
[Fell] luminosity at the optical nucleus is log[FeII] (erg s"-*^) < 33, far below the expected luminosity of a ~ lOMyr population. In the right 
panel, the region incompatible with the low Br7 and [Fell] fluxes is shown in grey. The blue strip shows the range of C02.22-2.37^m index 
allowed by our data to la confldence. 



6.1. Age and Mass Estimation from Stellar Kinematics 

We have extracted an integrated spectrum of the core 
of the nucleus using an aperture of radius i?eff ^ and per- 
formed a similar analysis as for the stellar kinematics, 
making use of the pPXF code to fit it with a library of 
stellar templates. The results of the fitting yield a veloc- 
ity dispersion of cr^^^ c± 71ibl4kms~^. We have adopted 
this value to provide an estimate of the dynamical mass 
of the cluster, given by the following relation: 



Mdyn T] 



G 



(1) 



where 77 is a geometric constant, <j is the velocity dis- 
persion, i^eff is the effective radius of the cluster and G is 
the gravitational constant. The geometrical constant 77 is 
determined by the density and velocity distribution of the 
cluste r, and is comnionly assumed to be ^10 when R ^^ is 
used (|Walcher et al.|[2005l : [McCrady & Graham"2007j to 
obtain the total mass of the cluster. However, different 
values from 3 to 10 co uld be adopted (see Hagele et al. 
I2QQ9L iBarth et al.|[2QQ9 and B asu-Zych et a l. 2009). As 
shown in Fig. El the optical nucleus has its own internal 
velocity gradient, distinct from the main velocity gradi- 
ent of the galaxy, that suggests that the cluster is viri- 
alized. Using the lowest value of 77 = 3 that assumes an 
isotropic velocity field, yields Mdyn = (1.1±0.4) x 10^ Mq 
for the mass of the cluster, a value consistent wit h an e ar- 
lier stellar kinematic estimate by .Thatte et al.l (|2000[ ) . 

Using both NICMOS F222M image and our SINFONI 
data, we have obtained the K band luminosity of the 
star cluster. In order to perform synthetic photome- 
try in our IFU data, we made use of the K band re - 
sponse curve of 2MASS as defined in lCohen et all (|2QQ1 
to obtain a synthetic K-band image. We have used 
the same aperture of radius 2.97±0.15pc to measure 
the luminosity of the star cluster in both images. The 
values obtained are i^K,siNFONi = 6.03 x 10^ and 
i^K,NiCMOS = 5.43 X 10^ Lq from the SINFONI and NIC- 
IVlOS data respectively. 

1 Lr^ = 3.826 X 10^6 W. 



We can thus derive a mass-to-light ratio Mdyn/L ra- 
tio based on the dynamical mass obtained above, which 
is Mdyn/L - 0.29 ± 0.12 Mq/Lk,q (where the denom- 
inator is in units of monochromatic solar K-band lumi- 
nosity I/K,0 = 2.150 X 10^^ W). In order to make a first 
estimation of the age of the cluster, we have us e d the 
stellar population synthesis models of 'Marastor] (|2005L 
MOb hereafter). The age that corresponds to the mass- 
to-light ratio derived above according to these models is 
logT(yr) 



-1 7+0.08 
' -0.36- 



6.2. Age and Mass Estimation From Spectral 
Diagnostics 

In addition to the age derived from stellar kinematics 
and photometry, we can put further constrains on the 
age via the CO stellar absorption. We have measured 
the equivalent width of the first CO bandhead as well as 
its CO index (defined as the ratio of the flux densities 
at 2.37/im and 2.22/im) to compare the values with the 
prediction from the SB99 and the 1VI05 synthesis models 
respectively. The value obtained for the CO equivalent 
width is Wco = 11.06 ± 0.25 A, where uncertainties are 
calculated by a IVIonte Carlo method of TV = 1000 simu- 
lations of the spectrum. 

To compare the equivalent width with the SB99 mod- 
els (assuming a Salpeter IIVIF and solar metallicity), we 
have considered an instantaneous burst of star forma- 
tion as implied by the low Br7 emission from the clus- 
ter. It is compatible both with a young population of 
less than logT(yr) ~ 7.17 or with an evolved population 
older than logT(yr) ~ 7.70. On the other hand, accord- 
ing to the models, the [Fell] emission suggests a lower 
limit of logT(yr) ~ 7.55 for the age of the star cluster. 
Combining the constraints from the [Fell] emission and 
the Wco yield an estimate for the age of the cluster of 
logr(yr) = 7.97i°;ii 

In contrast to SB99, the 1VI05 models use the CO in- 
dex. The integrated spectrum of the optical nucleus 
yields CO2.22-2.37/im = 0.082 ± 0.018. Using the M05 
models with the same constraints as discussed above 
from the Br7 and [Fell] emission, the resulting age is 
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logT(yr) = 7.98t°;^?. 

It is notable that the ages estimated here from the 
CO equivalent width and index are very similar to that 
calculated previously from the mass-to-light ratio using 
the dynamical mass. Instead of comparing ages derived 
using the different methods, we can instead compare 
masses. To do this, we first estimate the mass-to- light 
ratio M'^ / L from the M05 model associated with the age 
derived from the CO index. We can then derive the stel- 
lar (rather than dynamical) mass for the core of the star 
cluster. Again adopting a Salpeter IMF and solar metal- 
licity, we find W / L 0.22 ± 0.08. Since the luminos- 
ity is measured within i^eff (which by definition contains 
half the light), we set the luminosity of the cluster to 
be twice that given in Sec. 16.11 Thus the stellar mass is 
(7.8 ±2.4) X 10^ M0. 

This value is a factor of a few higher than the pho- 
tometric mass estimate from Thatte et al. (2000) and 
rather similar to that given in iWiklind et aD (r2004). We 
note that, within the inevitable uncertainties, it is also 
remarkably consistent with the dynamical mass we esti- 
mated previously. 

These ages (from 16.11 and 16. 2p are a little higher, but 
not inconsistent with that discussed at the end of Sec. [5] 
and all point to an age approaching, but perhaps a little 
less than, ~ 100 Myr. We have considered only instanta- 
neous star formation, since continuous models are clearly 
ruled out by the low Br7 equivalent width. However, a 
short but finite burst length of, for example, ~ 50 Myr 
could plausibly reconcile the small differences between 
the estimates. 

6.3. Could the Optical Nucleus Host a Supermassive 
Black Hole? 

As a particularly massive star cluster, and potentially 
the nucleus of M83, the optical nucleus is a suitable 
candidate for hosting a supermassive black hole. We 
have therefore estimated the mass of the supermassive 
black hole that one might expec t to find in the i nner 
regions of M83. Using GALFIT (|Peng et al.llMol ). we 
fitted a three component mo del to the 2MASS K-band 
image (Skrutskie et al. 2006), that allows us to sepa- 
rate the contribution of the disc, bar and bulge to the 
total flux distribution. The K-band luminosity of the 
bulge derived from the fit is logLx, bulge = 9.705 in Lkq 
units. Taking into ac c ount the Mbh/^k, bulge relation 
from Ma rconi fc HuntI (|2003[ ), we estimate the mass of 
the central black hole should be around ~ 3.9 x 10^ Mq. 
This value for the BH mass is al so similar to that ob - 
tained from the Mbh - cr relation (i Tremaine et al.|[2002[ ), 
if we take lOOkms"^ from Fig. [5] to be the velocity dis- 
persion of the stars in the bulge. Although this mass 
is less than the dynamical mass of the optical nucleus, 
it seems unlikely that a supermassive black hole should 
make up more than ~ 35% of the dynamical mass of a 
star cluster. 

On the other hand, the measurements from Chandra 
presented in Soria & Wu (2002) show that the optical 
nucleus is one of the brightest sources in X-rays within 
the nuclear regions of M83. They fit the nuclear spectra 
to a power-law model with total X-ray luminosity (0.3- 
8keV) of Lx ~ 2.6 x lO^^ergs"^ and a photon index 
of F ~1.15 (see their Table 3). The total luminosity is 



compatible with the X-ray emission from either a low 
luminosity AGN or a stellar-mass black hole candidate 
in a binary system (|Grimm et al.l l200"3l ). However, the 
lack of any other AGN signature (e.g. [SiVI] emission at 
1.96 /im) indicate that if a supermassive black hole were 
present, it would not be in an AGN phase at the present 
epoch. It is likely then that the X-ray emission comes 
from X-ray binaries in the hard state. 

7. LOCATION OF THE NUCLEUS 

The location of the nucleus of M83 has been topic of 
an int ense debate over the last decade. iThatte et al.l 
(|2000[ ) first reported the discovery of a double nucleus 
in M83 based on long-slit measurements of the velocity 
dispersion of stars. They found two peaks in the veloc- 
ity dispersion, one coincident with the optical nucleus, 
and another associated with the centre of symmetry of 
the bulge isophotes, both enclosing a dynamical mass of 
~ 1.3 X 10^ M0. However, the lack of two-dimensional 
information did not allow them to obtain the precise po- 
sit ion_oftln^ second mass concentration. 

Mast et al] (|2QQ6f ) made optical integral field spectro- 
scopic observations of the inner 12'' x 21'' of the galaxy 
that allow them to study in more detail the velocity field 
of the ionized gas. Their results support the conclusion 
about the presence of a second mass concentration of 
~ 1.0 X 10^ Mq (although it is possible that its position 
was misplac ed due to an incorrect spatial scaling) . 

Diaz et al. (2006) pointed out the position of a mass 
concentration at the northernmost part of the southern 
star-forming arc. They estimated that this hidden con- 
centration would enclose a total mass of ~ 1.6 x 10^ M©, 
derived from the ionized gas kinematics. They also found 
that the position of this hidden mass corresponds to a 
peak of emission in the mid-IR continuum at 10 jam. 

More recently, iHoughton fc Thatt^ (|2008[ ) combined 
near-IR long-slit spectroscopy with HST imaging to 
study the 20" x 20" central region of M83. They looked 
at the stellar kinematics for dynamical signatures of pu- 
tative hidden mass concentrations at locations indicated 
in previous work. Their results show no evidence of ob- 
scured masses and they conclude that the velocity gra- 
dients observed in the gas kinematics are a consequence 
of shocks. 

Making use of near- IR integ!:ral field spect roscopy and 
numerical simulations, Rodrigues et al. (2009) presented 
a detailed study of the ionised gas kinematics in the inner 
5"x 13" of the galaxy, covering the wavelength interval 
from 1.2 /im to 1.4 /im with a spatial sampling of 0'.'36. 
They focused on the dynamical properties and evolution 
of the optical nucleus, the CO kinematic centre (similar 
to the photometric centre) and the putative mass concen- 
tration proposed by iDiaz et al.i (2006) that is coincident 
with the 10/im continuum emission peak. Based on the 
ionised gas kinematics, they derived dynamical masses of 
- 6.0 X 10^ Mq, - 4 X 10^ Mq and - 2.0 x 10^ M© for 
the kinematic centre, the optical nucleus and the hidden 
mass concentr ation respe ctively. 

Finally, Kn apen et al.l jlOlO) analyzed in detail the 
ionized gas kinematics, making use of Pa/3 IFU obser- 
vations with a s patial sampling of Q '.'36, confirming the 
results found bv lThatte et al.l ()2000f ) that the photomet- 
ric centre coincides with the kinematic centre, and this 
location is offset by ~4" (~90pc) from the optical nu- 
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cleus. They proposed two possible options for the lo- 
cation of the true nucleus of M83. One option is the 
presence of an obscured hidden mass in the kinematic 
and photometric centre, that would require a dust ex- 
tinction of Ay = 3 — 10 mag. However, authors consider 
this option unlikely, because no other signatures - such 
as a peak in the velocity dispersion or in the near-IR 
emission - of a hidden mass are found. Their other op- 
tion was that the optical nucleus is the true nucleus, and 
it is displaced from the kinematic centre as a result of 
some past interaction. 

Our results clearly show that on scales of tens of par- 
sees, the gas kinematics are dominated by shocks and 
outflows. But we have also seen there is a global gradient 
of ~60kms~^ along the inner arc which is totally unre- 
lated to the stellar velocity field. Given the orientation 
of the galaxy, these gas kinematics appear to be tracing 
an infiow of gas along the inner bar to the photometric 
centre. We would expect such a bar-driven gas inflow 
to terminate at the nucleus, and the location of the end 
of the inflow does coincide with the photometric centre. 
We also emphasize that there is an increase of the stellar 
continuum at this point, consistent with the centre of the 
bulge. But this continuum, primarily from older stars, 
is mostly swamped by the stronger continuum from the 
young stars in the surrounding star-forming ring. The 
lack of a very recent starburst in the photometric centre 
would explain the low stellar continuum emission in this 
region. 

As mentioned before, we have not found any evidence 
or signature of a hidden mass near the kinematic centre. 
However, the velocity dispersion expected due to ran- 
dom motions around such a black hole, given the spatial 
resolution of our AO data, would be (Jbh ^ 60kms~^. 
Taking into account the instrumental resolution of our 
data {R ~ 1500, or an instrumental broadening of 
c^instr ~ 85kms~^) and the dispersion of the stars of 
cr^ ~ lOOkms"^ measured almost everywhere, the ef- 
fective enhancement of the dispersion expected to be ob- 
served for this BH would be of ~10kms~^, making our 
data insensitive to such a compact object. Thus, the lack 
of a kinematic signature of a supermassive BH does not 
equate to the absence of such a BH at the location of the 
kinematic and photometric centre of M83. 

We also considered the possibility that the ~ 3 A" off- 
set of optical nucleus from the kinematic centre could 
be the result of an m = 1 pertur bation in the ^ravita - 
tional potenti al, as sug:g;ested bv iKnapen et al.l (|2010l ). 
As shown in iBournaud et al.l ()200 5'). an m = 1 per- 
turbation of the potential could be explained by asym- 
metric accretion of gas towards the inner regions of 
the galaxy, which would be in agreement with the un- 
perturbed spiral pattern revealed at mid-infrared wave- 
lengths fo r M83 dPale et al.ll2QQ9D. A lternativelv. as dis- 
cussed in iHopkins fc OuataertI ()201Q| ) , an eccentric pat- 
tern resulting from a past interaction could persist up to 
~ 10^ dynamical times. However, a model of an eccentric 
disk as used by iTremaind (| 19951 ) to explain the double 
nucleus of M31 is probably not appropriate to explain the 
offset in here. It is caused by the high density of stars 
near apocenter in their elliptical orbits: the combination 
of their slow velocities at this point together with the fact 
that their velocity is along the line of sight, create the 



illusion of a secondary nucleus. In M83, the compactness 
of the optical nucleus and the presence of a coherent in- 
ternal velocity gradient (see Fig. [6]) are inconsistent with 
the rather diffuse appearance expected for such a grav- 
itational perturbation. Furthermore, the presence of an 
apparent secondary nuclei as in M31 is favoured by the 
close edge-on sightline. In principle, the eccentric disk 
could obscure the radiation from the black hole, even at 
X-ray wavelengths, and would explain the lack of signa- 
ture of a compact object near the kinematic centre of 
M83. However, the low inclination of the galaxy and the 
fact that the optical nucleus lies within the rather sym- 
met ric circumnuclear ring traced by the molecular gas 
(see iSakamoto et al.ll2004l i make it unlikely. 

8. CONCLUSIONS 

We have presented new and detailed near infrared 
adaptive optics integral fleld spectroscopic data for the 
innermost ~200pc of M83, and studied the kinematics 
and distribution of the stars as well as molecular and 
ionised gas. Our conclusions are as follows: 

1. The stellar kinematics show a smooth global veloc- 
ity fleld typical of uniform rotation. They reveal 
an independent and coherent velocity gradient in- 
trinsic to the optical nucleus with an amplitude 
of ~30kms~^. The velocity dispersion across the 
whole region also shows a smooth distribution, with 
values in the range ~ 50-100 km s~^. 

2. The ionized and molecular gas reveal a complex 
situation in which their kinematics are completely 
dominated by shocks and inflows at small scales but 
trace globally an inflow along the nuclear bar to the 
kinematic centre. This, and the fact that they are 
totally unrelated to the stellar kinematics, make 
them unsuitable to estimate dynamical properties 
of the central regions. 

3. There is plentiful evidence for recent super novae. A 
bright ring-like Br7 feature, which dominates the 
inner star- forming arc, has low dispersion and no 
measurable velocity gradient. It can be explained 
in terms of a light echo from a recent supernova 
explosion. The [Fell] emission has high dispersion 
along the arc suggestive of shocks from supernova 
remnants; and in the optical nucleus shows two 
bright locations where the dispersion is high, which 
are likely to be individual supernova remnants. A 
comparison of the gas and stellar kinematics in- 
dicates that the off-nuclear mass concentrations, 
which had been proposed on the basis of ionised 
gas kinematics, are instead regions where there are 
complex kinematics associated with recent super- 
nova events. 

4. A spatial study of the excitation mechanisms of the 
warm H2 suggest that the inner arc is dominated by 
collisional mechanisms, consistent with an episode 
of recent supernova events, while H2 in the optical 
nucleus has a higher contribution from radiative 
processes. 

5. The ~ 3pc effective radius of the optical nucleus 
together with its coherent internal kinematics yield 
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a dynamical mass of Mdyn = (1.1 ± 0.4) x lO'^ Mq. 
Its K-band luminosity is Lk = 5.7 x 10^ L©. 
The resulting mass-to-light ratio implies an age of 
100 Myr that is fully consistent with that im- 
plied independently by the CO index and equiv- 
alent width. Similarly, the age implied by the 
CO index yields a stellar mass for the cluster of 
W (7.8 zb 2.4) X 10^ M©, consistent with the 
dynamical mass. 

6. We show that the optical nucleus cannot be an 
m = 1 perturbation, and is not the 'true' nucleus of 
M83. Instead, we argue the this is indeed located 
at the photometric and kinematic centre of M83's 
bulge, where there is a measurable peak in the K- 
band continuum (albeit swamped by the bright sur- 
rounding star forming ring). We also show that, for 
the expected black hole mass in the centre of M83, 
one would not expect to see a dynamical signature 
at currently attainable spatial resolutions. 
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